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ABSTRACT 



Context. The dynamical evolution of a single globular cluster and also of the entire Galactic globular cluster system has been studied 
theoretically in detail. In particular, simulations show how the 'lost' stars are distributed in tidal tails emerging from the clusters. 
Aims. We investigate the distribution of Galactic globular cluster stars on the sky to identify such features like tidal tails. The Sloan 
Digital Sky Survey provides consistent photometry of a large part of the sky to study the projected two dimensional structure of the 
17 globular clusters in its survey area. 

Methods. We use a color-magnitude weighted counting algorithm to map (potential) cluster member stars on the sky. The same 
algorithm was already used by Odenkirchen et al. (2003) for the detection of the tidal tails of Palomar 5. 

Results. We recover the already known tidal tails of Pal 5 and NGC 5466. For NGC 4147 we have found a two arm morphology. 
Possible indications of tidal tails are also seen around NGC 5053 and NGC 7078, supporting earlier suggestions. Moreover, we find 
potential tails around NGC 5904 and Pal 14. Especially for the Palomar clusters other than Pal 5, deeper data are needed in order to 
confirm or to rule out the existence of tails. For many of the remaining clusters in our sample we observe a pronounced extratidal 
halo, which is particularly large for NGC 7006 and Pal 1. In some cases, the extratidal halos may be associated with the stream of the 
Sagittarius dwarf spheroidal galaxy (e.g.,NGC 4147, NGC 5024, NGC 5053). 

Key words, globular clusters: general - globular clusters: individual: NGC 2419, NGC 4147, NGC 5024 (M53), NGC 5053, 
NGC 5272, NGC 5466, NGC 5904 (M5), NGC 6205 (M13), NGC 6341 (M92), NGC 7006, NGC 7078 (M15), NGC 7089 (M2), 
Palomar 1, Palomar 3, Palomar 4, Palomar 5, Palomar 14 
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1. Introduction 

Globular clusters (GCs) are the oldest stellar objects found in the 
Milky Way (MW), with ages in the range of 10 - 13 Gyr. They 
represent tracers of the early formation history of our Galaxy. 
GCs are ideal systems to study stellar dynamics in high-density 
systems as their relaxatio n times are much smaller than their age 
dGnedin & Ostrikerf 19971 G097). We expect, at least in the clus- 
ter core, the stars to have lost memory about their initial condi- 
tions. The GCs we see today are only the lucky survivors of an 
initially much larger population (G097). Several different mech- 
anism s can dissolve a GC: internal processe s include infant mor- 
tality dGever & Burkerd200lHKroupal200Th by stellar evolution, 
gas expulsion and two-body relaxation; external mechanisms in- 
clude primarily the interaction of the GC with the constantly 
changing tidal field of its host galaxy. Different processes are 
dominant at different ages of the cluster. For young clusters stel- 
lar evolution is the main driver that may potentially cause them 
to dissolve. More massive stars are evolving rapidly, have strong 
winds and explode as supernovae within the first few million 
years driving out the remaining gas. The sudden change, due to 
the mass loss, in the GCs internal gravitational field can destroy 
the GC. Two-body relaxation arises from close encounters be- 
tween cluster stars, leading to a slow diffusion of stars over the 
tidal boundary. This process is taking place from the beginning 
to the final dissolution of the cluster. Equipartition of energy by 
two-body encounters leads to mass-segregation, the more mas- 



sive stars sink to the center, less massive stars remain in the outer 
parts of the cluster. This leads to a preferential loss of low-mass 
stars, which ca n be observed today in decl ining mass-function 
slopes of GCs (Baumga rdt & Makinol l2003). On the other side, 
the mass loss due to the cluster's encounter with the Galactic 
disk and bulge is stronge st at pericenter passages, esp ecially for 
GCs on elliptical orbits (Baumgardt & Makino 2003). All these 
mechanisms have one common result: the GC is losing stars. As 
a result of relaxation stars gradually diffuse out, while during 
tidal shocks the GCs loss of stars is more sudden. 

1.1. Observational studies 



iGrillmair et alj d!995l) performed deep two-color star counts 
to examine the outer structure of 12 GCs using photographic 
Schmidt plates. They detected a halo of extratidal stars around 
most of their clusters. They also found indications of possible 
tidal tails emerging from some clusters. The tidal tails are not in 
all cases purely of stellar origin, as the morphological identifi- 
cation of galaxies and stars, especially for faint objects, was not 
perfect. I.e., in some cases background galaxies were mistaken 
for stars. In other cases random overdensities of foreground stars 
induced features into the cluster's contours. [Leon et aT I d2Q00T > 
studied the tidal tails of 20 GCs, using Schmidt plates as well. 
These authors also find halos of extratidal stars as well as tidal 
features for most of their clusters. The two stu dies have three 
GCs in common. For two out of the three clusters Grillmair et al. 
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(1 19951) detected tidal tails where iLeon et al l (l2000h did not see 
such a strong signal and in the third case it is the opposite. 

Apart from the two studies mentioned above, different au- 
thors determined the 2d distribution of cluster stars on the sky 
for individual clusters. The two most prominent cases are Pal 5 
and NGC 5466, as for these two clusters extended tidal tails have 
been found. The large area-coverage with CCDs by the SDSS 
provides a large database for tidal tail studies of GCs. Pal 5 was 
the first cluster for whi ch extended tidal t ails were discovered 
and s t udied in detail (jO denkirchen et al. 2001; iRockosi et alj 



2002t lOdenkirchenetaiT 
Odenkirchen et all l2QQ9h 
was the second cluster 



2003 



tails 



iGrillmair & Dionatosl 120061: 
using SDSS data. NGC 5466 
fo und t o have extended tidal 
(lOdenkirchen & Grebell l2004t iBelokurov etail 120061: 



IGrillm air & Johnson! 120061) . Not only spectacular tidal tails, 
but simpler tidal feat ures were also detected for NGC 5053 
dLauchner et alJ 120061) . for N GC 6341 (M92) bv iTesta et al.1 
(2000|), and for NGC 414 7 by iMartfnez Delgado et al.l (120041) . 
Further, iKiss et aL I J2001 studied the kinematics of red gi- 
ants around five GCs. For M55, M30 and NGC 28 8, they did 
not fin d any strong signs of tidal debris. Recently, IChun et all 
(2010) studied the spatial configuration of stars around five 
metal-poor GCs (M15, M30, NGC 5024 (M53), NGC 5053, and 
NGC 5466). They used deep images obtained at the Canada- 
France-Hawaii Telescope. Around all GCs extratidal overdensi- 
ties and extratidal halos were observed. Between NGC 5024 and 
NGC 5053 they detect a tidal-bridge like feature and an envelope 
structure. 
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Fig. 1. Vital diagram from G097. The red stars are the GCs in 
our sample, Cluster mass and half-light radius were taken from 
G097. GCs within the triangles are most likely surviving the 
next Hubble time. Adopted from G097, reproduced by permis- 
sion of the AAS 



1.2. Theoretical prediction and N-body simulations 

iGnedin & Ostrikeri dl997l) (G097) studied the dynamical evo- 
lution of the entire Galactic GC system. They derived the de- 
struction rates due to the internal and external mechanisms. 
According to their calculations, the total destruction rate of the 
entire GC system is such that more than half of the present day 
GCs will not survive the next Hubble time. G097 did not in- 
clude any information on the GCs' proper motion, they used the 
known present day distance to the Galactic center and the GCs' 
rad ial velocity. 

iDinescu et"aTI (1 1999b (DG99) derived the destruction rates 
for a smaller sample of Galactic GCs in a similar way, but in- 
cluded the measured proper motion data. They concluded that 
the destruction processes for the clusters in their sample are 
mostly dominated by internal relaxation and stellar evaporation. 
Tidal shocks due to the bulge and disk are only dominant for 
a small number of GCs in their sample. For the sample we are 
studying here this is only the case for Pal 5. Comparing their 
destruction rates with the numbers of G097, they find that the 
G097 total rates are larger. The destruction rates due to internal 
processes (2-body relaxation, stellar evaporation) are compara- 
ble in both studies, therefore only the destruction rate due to 
tidal shocks is incompatible. DG99 conclude that the statistical 
approach of assigning tangential velocities (the method used by 
G097) results in more destructive orbits than are actually ob- 
served. 

A similar conclusion is drawn by lAllen et aL (2006), who 
derive the orbits of Galactic GCs, also based on observed proper 
motions, in an axisymmetrical potential and in a barred poten- 
tial (resembling the MW). They derived the destruction rates for 
the GCs in their sample in an almost identical way to DG99. 
lAllen et ail (120061) further investigated the influence of the bar 
on the clusters' orbits and destruction rates. The bar does not 
influence the orbits of clusters that have a pericenter distance 



greater than ~ 4 kpc. Hence, none of the clusters in our sample 
(see Table [TJ are influenced by the bar. 

ICombes et al.l (1 19991) studied the tidal effects experienced by 
GCs on their orbits around the MW. In general, they find that 
two large tidal tails should emerge extending out to 5 tidal radii. 
They conclude by stating ;) GCs are always surrounded by tidal 
tails and tidal debris, ii) the tails are preferentially composed of 
low mass stars, Hi) mass loss in a GC is enhanced if the cluster is 
rotating in the direct sense with respect to its orbit, iv) extended 
tidal tails trace the cluster's orbit, and v) stars are not distributed 
homog eneously along the tidal tails, but clump. 

Canu zzo Dolcetta et alJ (120051) investigated the clumpy 
structure of the tidal tails. They found that the clumpy structure 
is not associated with an episodic mass loss or tidal shocks due to 
encounters with Galactic substructure, as stars are continuously 
lost. T he clumps are not self-gravitating systems. iKiipper et al.l 
(2008) showed that the clumps ar e a result of ep icyclic motion 
of the stars lost by the cluster. lJust et all (120091) investigated 
the formation of th e clumps and discussed their observability. 
Kiip per et alJ d2009t) investigate the position of the clumps rela- 
tive to the cluster for clusters on eccentric orbits. They concluded 
that tidal shocks are not the d irect origin of the clumps. Further, 
Capu zzo Dolcetta et al.1 d2005l) confirmed the earlier finding that 
the cluster's leading tail develops inside the orbit, while the trail- 
ing tail follows outside the orbital path. They also found that for 
clusters on elliptical orbits the tid al tails trace the or bital path 
only near perigalacticon (see also iDehnen et alJl2004l) . On the 
other hand for circular orbits the tails are clear tracers of the 
cluster path. Moreover the length of the tails is not constantly 
growing. On eccentric orbits the leading tail tends to be longer 
than the trailing tail during the cluster's motion from apogalacti- 
con to perigalacticon and vice versa during the other half of the 
cluster's orbit. 
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iMontuori et al.l (120071) investigated the direction of tidal tails 
with respect to the GC's orbit. In the outer parts of the tails (> 
7-8 tidal radii away from the cluster center), the tails are very 
well aligned with the cluster's orbit regardless of the cluster's 
location on the orbit. On the other hand, in the inner parts, the 
orientation of the tidal tails is strongly correlated with the orbital 
eccentricity and the GC's location on the orbit. Only if the cluster 
is near perigalacticon, the inner tidal tails are aligned with the 
orbital path. Therefore, only if long tidal tails are detected, it 
is possible to securely constrain the cluster's orbit from those. 
Detecting only small, short tidal extensions just outside the GC's 
tidal radius do not give any hint on the cluster's orbit unless the 
cluster's proper motion has been measured before. 

DG99 published proper motions of a large number of GCs 
and added values from the literature to their catalog. From this 
they derived orbital parameters, such as ellipticity e, perigalacti- 
con Rp er i, apogalacticon R apo , etc. The eleven clusters we have in 
common with DG99 have orbital eccentricities larger than ~ 0.3, 
i.e., there are no circular orbits for our sample. 

These studies give a good theoretical understanding of the 
size, shape, orientation, etc. of tidal tails of GCs. In an observa- 
tional stud^^fjigji^£lement_ abundance variations in field halo 
stars, Martell & Grebel (2010) conclude that as many as half 
of the stars in the Galactic halo may originally have formed in 
meanwhile dissolved GCs. In this paper we use data from the 
Sloan Digital Sky Survey to study the 2d structure of 17 Milky 
Way GCs and to compare the observed features with the theo- 
retical predictions. The paper is organized as follows: Section[2] 
introduces the photometric data used in our analysis. Section [3] 
explains the algorithm used to count (potential) member stars on 
the sky. In Section [4] we present the number density profiles of 
the 17 clusters in our sample. Section[5]presents the derived con- 
tour maps for our GCs. In Section[6]we discuss our results. The 
paper is concluded in Section|7]with a summary. 



merged the two datasets into one catalog for each GCs. The mag- 
nitudes in the final catalog were corrected for extinction. The 
SPS S catalog provides the Galactic foreground extinction values 
from lSchlegelet all dl998h . For the fields with An08 photome- 
try, the extinction values were derived by a cubic interpolation 
of the values listed in the SDSS catalog. 

In TableQ]we list the GCs found in the SDSS DR7 footprint. 
Columns (3) and (4) contain the equatoria l coordinates ( J2000.0) 
of the clust er center mai nly taken from lHarrii d!996|). except 
for P al 14 dHilkerf [2006h and NGC 7089 dDalessandro et al. I 
2009). In column (5) we list the clusters' distance to the Sun 
and in col umn (6) the c lusters' distance to the center of th e 
MW from lHarrisI dl996h . except for Pal 14 dHilkerl [2006). 
Columns (7) and (8) contain the core and tidal radii from 
McLaughlin & van der Marel ] d2005h (M05). For almost all GCs 
in our sample proper motions have been measured. Only for 
three of the most distant GCs, NGC 2419, Pal 4 and Pal 14, as 
well as for NGC 5053 this information is lacking. Columns (9) 
and (10) list the proper motion taken from DG99. 

In Figure [2] we show the color-magnitude diagrams (g-r,g) 
of all GCs in our sample in order of increasing distance from 
the Sun (see column 5 in Table [TJ. We show all stars within the 
clusters' tidal radii. In the bottom panels we show for each GC 
a sample of field stars. The field stars are chosen from a random 
field in the vicinity of the GC (~ 2° away) within a circle equal 
to the clusters' tidal circle. 

The SDSS data were taken over several years and in 
differ ent observing conditions. The average seeing is 1.5" 
(lAdelman-McCarthv et al. 20p7|). In our study we only use the 
g, r, and i band. Ivezic et al. I ( 20041) showed that stellar popula- 
tion colors are constant over the survey area, concluding that 
the photo metric calibration is accurate to roughly 0.02 mag. 
Further, llvezic etaLl (120071) analyzed multiple observations of 
bright stars in the southern equatorial stripe. The photometry was 
repeatable to less than 0.01 mag in all five bands. 



2. Data 

The Sloan Digital Sky Survey (SDSS) is an i maging and 
spect roscopic survey in the Northern hemisphere dYork et alj 
2000). SDSS 



imaging 



ugriz (Fukugitaet al 
200 1: ISmith et alj|2002tlivezic et 



data a re produced 
1996; iGunn 



produced in th e five 1 
inn et al.l [1991 iHogg 
al. 1 12004 lfuckeret al.l l 



five bands 
etal 

200< 



Gunn et al.ll2006l) . The data are automaticall y processed to mea- 



sure p hotometric and astrometric properties (lLupton et alj 2002; 
iPier et al. l[2003l Photo) and are publicly available on the SDSS 



web pagefl We used the da ta from the SDSS Data Release 7 
(DR7; lAbazaiian et al.ll2009l) . 

The automatic SDSS pipeline Photo was initially designed 
to process high Galactic latitude fields with a low density of 
Galactic field stars. Fields centered on and around GCs are too 
crowded for Photo to process, so the automatic pipeline does not 
provid e photometry for these most crowded regions. lAn et al.l 
(20081 An08) used the DAOPHOT crowded-field photometry 
package to derive accurate photometry for the stars in these 
crowded areas. This photometry is published on the SDSS web 
pages as a value-added catalogs From this catalog we only con- 
sidered stars with photometry flags set to 0. The An08 catalog 
and the original SDSS catalog overlap around each GC. If a star 
is measured in both samples, we used the original SDSS pho- 
tometry. This choice has no influence on the result. Finally, we 



1 www.sdss.org 

2 http://www.sdss.org/dr7/products/value_added/index.html 



3. Counting Algorithm 

To map c luster member star s on t he sky we used the same 
method as lOdenkirchen et"aL (|2003[) to follow t he tid al tails of 
Pal 5. We refer the reader to lOdenkirchen et al] d2003l) for a de- 
tailed description of the procedure. Here, we only give a rough 
overview. 

For the subsequent study of the distribution of stars on the 
sky we need a clean sample of potential member stars. In order 
to reduce the contamination of our sample with field stars we 
define new orthogonal color indices: 



c\ = k l -(g-r)+k 2 -(g- i) + h 
c 2 = -k 2 ■ (g - r) + ki ■ (g - i) + k A . 



(1) 

(2) 



In a color-color plot the cluster stars are located in an almost 
one-dimensional locus. The new indices were now chosen in 
such a way that the color c\ traces the cluster stars' locus and c 2 
is perpendicular to that. To derive these indices for each GC we 
used the stars within 2/3 • r,, where r, is the cluster's tidal radius 
from Table Q] We fitted a straight line y = tan(o') ■ x + b to the 
distribution of stars in the color-color diagram (g - r) vs. (g - /), 
where x represents the (g-r) axis and y the (g—i) axis. In Table|2] 
we list the fitted parameters tan(a) and b from which we derived 
k\ - cos(a), k 2 = sin(o-), ki — -b ■ sin(a), and k\ — -b ■ cos(a) 
for each GC. 

We selected our stars in the color-magnitude diagrams (c\ , g) 
and (c 2 , g) to minimize the contamination by Galactic field stars. 
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Table 1. Globular clusters in the SDSS DR7 footprint. 



NGC 


Name 


RA 


Dec 


"0 


□ 

"Af W 


r c 


ft 


(i a cos(<5) 








(hh mm ss.s) 


(° ' ") 


(knc) 


(knc) 




C) 


( mas/vr^ 


( mas/vr) 


2419 




07 38 08.5 


+38 52 55 


84.2 


91.5 


0.32 


8.34 






4147 




12 10 06.2 


+ 18 3231 


19.3 


21.3 


0.09 


6.32 


-1.85+0.82 


-1.3 + 0.82 


5024 


M53 


13 12 55.3 


+ 18 10 09 


17.8 


18.3 


0.35 


16.06 


+0.5 ± 1 


-0.1 ± 1 


5053 




13 16 27.0 


+ 17 41 53 


16.4 


16.9 


1.91 


13.85 






5272 


M3 


13 42 11.2 


+28 22 32 


10.4 


12.2 


0.37 


30.15 


-1.1 ±0.51 


-2.3 ± 0.54 


5466 




14 05 27.3 


+28 32 04 


15.9 


16.2 


1.53 


13.64 


-4.65 ± 0.82 


0.80 ± 0.82 


5904 


M5 


15 18 33.8 


+02 04 58 


7.5 


6.2 


0.42 


22.45 


+5.07 ± 0.68 


-10.7 + 0.56 


6205 


M 13 


16 41 41.5 


+36 27 37 


7.7 


8.7 


0.55 


19.49 


-0.9 + 0.71 


+5.5 ± 1.12 


6341 


M92 


17 17 07.3 


+43 08 11 


8.2 


9.6 


0.26 


12.96 


-3.3 + 0.55 


-0.33 ± 0.70 


7006 




21 01 29.5 


+ 16 11 15 


41.5 


38.8 


0.32 


12.75 


-0.96 ± 0.35 


-1.14 + 0.40 


7078 


M 15 


21 29 58.3 


+ 12 1001 


10.3 


10.4 


0.07 


21.50 


-0.95 + 0.51 


-5.63 + 0.5 


7089 


M2 


21 33 27.3 


-00 49 23 


11.5 


10.4 


0.32 


12.75 


+5.9 + 0.56 


-4.95 ± 0.86 




Pal 1 


03 33 23.0 


+79 34 50 


10.9 


17.0 


2.28 


12.04 








Pal 3 


10 05 31.4 


+00 04 17 


92.7 


95.9 


0.35 


4.89 


0.33 + 0.23 


0.30 + 0.31 




Pal 4 


11 29 16.8 


+28 58 25 


109.2 


111.8 


0.33 


3.30 








Pal 5 


15 16 05.3 


-00 06 41 


23.2 


18.6 


2.28 


12.04 


-1.78 + 0.17 


-2.32 ± 0.23 




Pal 14 


16 11 00.6 


+ 14 57 28 


74.7 


69.0 


0.77 


7.20 







Table 2. Parameters to derive the new color indices. 



NGC 


2419 


4147 


5024 


5053 


5272 


5466 


5904 


6205 


6341 


tan(tt) 
b 


1.472 
-0.039 


1.462 
0.004 


1.495 
-0.053 


1.434 
-0.006 


1.473 
-0.053 


1.531 
-0.047 


1.447 
-0.016 


1.381 
-0.010 


1.345 
-0.001 


NGC/Name 


7006 


7078 


7089 


Pal 1 


Pal 3 


Pal 4 


Pal 5 


Pal 14 




tan(cf) 
b 


1.540 
-0.028 


1.522 
-0.035 


1.473 
-0.054 


1.545 
0.007 


1.404 
0.017 


1.534 
-0.147 


1.553 
-0.034 


1.767 
-0.190 






Fig. 3. Color magnitude diagrams (ci,g) and (C2,g). Plotted are 
stars of NGC 6341 within one third of the cluster's tidal radius. 
The gray solid lines denote the borders within which the stars 
were selected. 



Figure [3] shows the CMDs (c u g) and (c 2 ,g) of NGC 6341. In 
(c2,g) we selected all stars with \c 2 \ < 2.5 ■ cr(c 2 )(g). <r(c 2 )(g) 



is the rms dispersion in c 2 for stars of magnitude g. Stars with 
deviating c 2 colors are unlikely to be cluster members. 

In (ci , g) we selected all stars within 3cr of the cluster's ridge 
line. The ridge line of a GC traces the location of highest density 
of stars in a color-magnitude diagram (CMD). Figure [3] shows 
the selection criteria for NGC 6341 as an example. This prese- 
lection leaves us with the final catalog, where all stars are explic- 
itly chosen to have the same photometric properties as the GC, 
but are located not only at the position of the cluster but also 
in the wide field around it. The final catalog is then split into the 
cluster sample to contain all stars within 2/3-r, and the field sam- 
ple to contain all stars outside 3/2 • r,. If a cluster was observed 
to be more extended than previously thought, we redefined the 
two samples. In some cases, e.g., for NGC 5024 and NGC 5053, 
two clusters are close neighbors in the projection on the sky. In 
these cases, we excluded all stars within the tidal radius of the 
neighboring cluster from the field sample. 

The stars in the preselected samples are not simply counted. 
The ratio of the distribution of cluster stars to field stars in the 
CMD (c; , g) was used as a weight in the sum. Figure|4]shows the 
two Hess diagrams of NGC 5272 as contour plots. The left panel 
is the field sample, the right panel the cluster sample. The cluster 
sample shows a prominent main-sequence and subgiant branch. 
Furthermore, a red giant and horizontal branch are visible. For 
the field sample, the density is more homogeneous in the entire 
selected areas. The horizontal branch is only prominent at the 
red edge. 

The density of cluster stars rtc(k) at a given point k on the 
sky is thus given as 
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Fig. 2. Top rows: C olor-magnitude diagrams of all GCs in our sample in order of increasing distance from the Sun. The distances 
from iHarrisI d 19961) are printed in each CMD. The grey horizontal line marks the magnitude limit for each cluster. Bottom rows: 
Color-Magnitude diagrams of a random field in the vicinity of the cluster (see text for details). 



n c (k) 



E; n (J, k) ■ PcU)/PfU) - n F (k) 

YjjP 2 C (j)/PF(j) 



(3) 



where pc and pf are the density distribution of the cluster and 
field sample, respectively, in the CMD (see Figure©, rif is the 
average background in the field around the GC. 

From Equation [3] we see that each star in a given cell on the 
sky k and in color-magnitude space j is weighted according to 
its position in the CMD by the factor pc(j)l 'Pf(J)- This weighted 
number of stars per cell on the sky is then summed and divided 
by the factor a = Y,j PcU)/PfU)- This gives the estimated num- 
ber of cluster stars nc plus a term rif/a, the number of contami- 
nating field stars attenuated by a. 

4. Number density profiles 

The large area coverage of the SDSS data allows us to study 
the number density profiles of GCs further out than any survey 
before. By using a large area "far away" from the cluster we are 
able to determine the mean background contamination and to 
investigate in detail the outer boundaries of the GCs. 

We derived the number density profiles by counting the stars 
in concentric annuli of constant width dr. We used only the stars 
preselected in colors and magnitude as described above. We used 



annuli of width 2' outside the tidal radius. Within the cluster's 
tidal radius we chose either 0.5' or 1' as width of the annuli, 
depending on the surface density of the cluster. An alternative 
way to derive the number density profile would be by using an- 
nuli containing a constant number of stars. Comparing these two 
methods revealed no differences in the resulting profiles and pa- 
rameters. 

For the final number density profile we counted all stars in 
each annulus and divided by the area of the annulus. The error 
of the number density was derived by dividing the annulus in 20 
segments. We then derived the number of stars in each segment 
and used the standard deviation of these 20 values as the error. 
For clusters with a small number of stars we divided the annuli 
in 8 segments to reduce the effect of small number statistics. 

For the denser clusters the SDSS is not suited to observe 
the number density profile in the cluster's central part due to 
crowding. Therefore, we combined the number density pro- 
files from the SDSS with previously published profiles for 
the inner pa r ts of t he GCs. The surface brightness profiles by 
Trag eretall dl995|) (T95) and the revised study of the same 
data by McLaughlin & van der Ma rel (2005) (M05) are the only 
large catalogs of surface brightness profiles of Galactic GCs. 
We combined these datasets with our new SDSS profiles to 
get profiles ranging from the cluster center out to several times 
the tidal radius. In order to combine the two datasets, we con- 
verted the surface brightness p(r) of T95 into a number den- 
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Fig. 4. Hess diagrams of the field (left panel) and cluster (right 
panel) sample of NGC 5272. Only stars already selected in Ci 
and C2 were used. 

sity log(N(r)) = -fi(r)/2.5 + C. C is a constant derived sepa- 
rately for each GC in the following way. The T95 profile and the 
SDSS profile overlap within a certain radial range. This over- 
lap was used to derive the vertical shift C between the profiles. 
The T95 profiles were shifted by C to match the SDSS profiles. 
For GCs whose T95 profiles extend far out, the outermost points 
were not taken into account, because these data points are the 
most affected by the background contamination or subtraction. 
Furthermore, differences in the outer parts of the different pro- 
files can also be explained by mass segregation in the clusters, 
unveiled by different limiting magnitudes. 

For the GCs NGC 5272 (M3), NGC 5904 (M5), NGC 6205 
(M13) NGC 6341 (M92) NGC 7 078 (M15), and NGC 7089 
(M2) iNovola & Gebhardtl d2006l) (N06) published surface 
brightness profiles derived from HST observations. They tested 
the influence of the chosen filter on the profiles and found no de- 
pendence. The profiles looked the same in each observed filter. 
For 50% of their clusters the photometric data points are brighter 
in the centra l parts of the clusters than the fitted King profiles 
(King][l962|). The discrepancy between the observations and the 
fit increases as the radius decreases. N06 conclude that the inner 
parts are not suited to be fitted with a King profile. 

For these six GCs we combined the N06 profiles and our 
SDSS profiles to span the largest possible range in radius. The 
combination was done in an identical way as for the combination 
of the T95 sample with our SDSS profiles. 

The c ombined profiles were fitted with iKingl (119621) and 
iPlummerl d 1 9 1 lb profiles. For the King profiles, we adopted the 
background levels derived from the SDSS observations and de- 
termined for the core and tidal radius through our fit. For the 
Plummer profiles we derived the Plummer radius r p through our 
fit. In Figures [5]4Z] we show the combined number density pro- 
files of all clusters in our sample. Table [3] lists the fitted core 
and tidal radii, the observed stellar background density tibkg and 
the Plummer radius. For some GCs, e.g. NGC5466, NGC 6205, 
the Plummer profile provides a better fit than the King profile. 
But in most cases, the Plumme r profile does no t trace the ob- 
served number density profile. Penarrubia et alJ d2009l) found 



that a Plummer profile more accurately traces the profile of a 
stellar system (mainly a dwarf spheroidal galaxy) that loses mass 
due to tidal effects. For NGC 5466 this is certainly the case and 
the known tidal tails are proof of tidal mass loss. For Pal 5, which 
is also known for its long tidal tails, the Plummer profile is not a 
better fit than the King profile. In the following Section we will 
discuss the prominent features and trends in the density profiles 
for each cluster individually. 

In Figures [814241 we plot our newl y derived t idal radii as red 
ellipses. The tidal radius, r t , from the lKingldl962l) profiles is the 
radius at which the stellar density (number of stars per arcmin 2 ) 
drops to zero. It must not necessarily be the radius outside of 
which a star has left the cluster potential, but it is an estimate 
of this radius. The radius outside of which the star left the clus- 
ter potential is the distance between th e cluster center and th e 
first Lagrangian point, the Jacobi radius ( Baumg ardt et alj2 009). 
Stars outside the Jacob i radius have left the GCs potential. 
Baumgar dt et alJ ([2009) calculated Jacobi radii for all Galactic 
GCs. We adopted these values. In Figures [814241 the Jacobi radii 
are plotted as blue dashed ellipses. 

5. Tidal Features 

In all contour maps we plot iso-density lines at levels of 
l,2,3,5,7,8,10,20,40,60,80,100,200,400,600,800,1000,2000-cr M? 
above the mean weighted background level n^ k . The mean 
weighted background was derived at least two tidal radii away 
from the cluster and if there were unusable areas, e.g., close 
neighboring clusters, overlapping SDSS scans, areas outside 
the SDSS footprint, known tidal tails, etc., these were not taken 
into account for its determination. <x Mg denotes the standard 
deviation of the weighted density in the area, which was used 
to determine n^ kg . The mean weighted background, n^ kg , is 
different from the background «m s derived in Section [4] n^g is 
the mean of the number of stars per area, whereas for the 
stars were weighted by the factor pc(j)l PfU) according to their 
position in the CMD. 

In Fi gures [514241 we present the contour maps for all 17 GCs 
in our sample. The upper/left panel shows a large area, typically 
9 deg by 9 deg, in order to show potential large scale features 
in a Lambert conical projection. The lower/right panels show a 
blow-up in order to study the small scale tidal structures of the 
GCs. On the x-axis of the zoom-in plots we plot the coordinate 
RAshown = X m +(RA-X m )-cos(Y m ), where (X m ,Y m ) is the center of 
the cluster. All contour plots have North up and East right. The 
calculated iso-density contours are smoothed with a Gaussian 
kernel over 5 pixels (= 15 arcmin) with a sigma of 2 pixels (= 
6 arcmin). In all maps the l<x Mg and the 2cr M? are shown in grey, 
while the higher contours are shown in black. 

Before presenting all contour maps, we briefly discuss the 
possible sources for contamination. 

- Background galaxies & quasars 

The SDSS automatic Star/Galaxy separation seems to be 
very reliable. As we only use data considered to be a point 
source we, in principle, do not have to worry about contam- 
inating background galaxies or galaxy clusters. In the fields 
with An08 photometry this may be an issue, but these fields 
do have a high stellar density, with stars outnumbering con- 
taminating galaxies. 

As quasars are point sources and are included in the 
SDSS database together with stars, these might influence 
our data. We have looked at the SDSS quasar catalog IV 
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Table 3. Measured structural parameters. 



JNLrL 


Name 




r c 


r t 








stars/arcmin 2 


arcmin 


arcmin 


arcmin 


2419 




0.35 


0.30 ± 0.05 


8.07 ± 0.87 


0.59 ± 0.01 


4147 




0.20 


0.14 ±0.01 


6.26 ± 1.53 


0.12 ±0.01 


5024 


M53 


0.22 


0.37 ±0.01 


14.79 ±7.19 


0.66 ± 0.01 


5053 




0.08 


2.23 ± 0.36 


12.52 ±4.49 


3.00 ± 0.3 


5272 


M3 


0.17 


0.41 ±0.003 


31.81 ± 12.3 


0.77 ±0.01 


5466 




0.03 


1.19 ±0.09 


14.63 ± 2.40 


2.32 ± 0.08 


5904 


M5 


0.52 


1.04 ±0.06 


18.91 ±2.71 


2.18 ±0.05 


6205 


M 13 


0.33 


0.82 ± 0.005 


14.36 ± 0.96 


1.27 ±0.01 


6341 


M92 


0.15 


0.37 ±0.01 


12.55 ±0.90 


0.36 ±0.01 


7006 




1.14 


0.14 ±0.01 


3.64 ± 0.56 


0.29 ±0.01 


7078 


M 15 


0.16 


0.27 ± 0.03 


18.35 ± 1.97 


0.26 ± 0.01 


7089 


M2 


0.13 


0.29 ± 0.002 


11.72 ±0.63 


0.35 ±0.01 




Pal 1 


0.13 


0.005 ± 1.06 


4.09 ± 0.92 


0.49 ± 0.03 




Pal 3 


0.03 


0.37 ± 0.02 


5.75 ±0.92 


0.65 ±0.01 




Pal 4 


0.15 


0.26 ±0.10 


5.30 ± 0.65 


0.54 ±0.01 




Pal 5 


0.15 


2.34 ±0.14 


16.03 ± 2.24 


3.38 ±0.11 




Pal 14 


0.07 


0.72 ± 0.08 


8.82 ± 3.20 


1.19 ±0.06 



dSchneider et al.ll2007l) . which contains 77 429 quasars from 
SDSS DR5. We applied the same selection criteria in color 
and magnitude to select possible contaminating quasars for 
each GCs, respectively. E.g., for NGC 5272 we end up with 
54 quasars contaminating our sample of 51 964 stars in total. 
The contamination due to quasars can be neglected. 

- Dust 

Extinction by dust might also influence our tidal features. 
All our magnitudes are corrected fo r extinction as given i n 
the SDSS database by values from ISchlegel et a"lf(ll998l) . 
In a first look, we do not expect a strong contamination by 
dust, as the SDSS mainly observed the sky around the North 
Galactic pole (NGP), where less dust is expected. But, in pre- 
vious studies a correlation between tidal tails and dust was 
found for some GCs. We examined for each GC the orienta- 
tion of the tidal features and studied the potential correlation 
with the dust. We did not observe any correlation. 

- Foreground & background stars 

In Section[3]we explained our criteria to select stars in color- 
magnitude space to minimize the number of stellar fore- 
ground and background contaminants. For all GCs we have 
only selected stars within 3<x of the cluster's ridge line to 
assure that the majority of stars counted are likely cluster 
members. The 2d-distribution of these selected stars shows 
in all cases a unique maximum at the GCs position. The area 
around the GC is in all cases not perfectly flat, but shows 
density peaks of various significances. I.e., we are able to 
clean our stellar sample, but not perfectly. A more accurate 
exclusion of Galactic foreground and background contami- 
nants can only be achieved with spectroscopic confirmations 
for each single star - a task which seems almost impossible 
for the time being. 

5.1. NGC 2419 

NGC 2419 is a very bright, unusually large GC located at a 
very large distance from the center of the Milky Way (Rmw = 
91.5 kpc, H96). Previous studies of its surfac e brightness pro- 
file revealed some evidence for extratidal stars (lBellazzinl 2007) 
and references therein). In Figure [5]our combined number den- 
sity profile is shown. We also observe a break in the profile's 
slope. No proper motion information is available. As NGC 2419 



is very remote, only the upper part of the red-giant branch and 
the horizontal branch are visible in the SDSS data. 

A stripe north of NGC 2419 has an artificially doubled 
background density. We corrected for this. All stars in the area 
denoted by the dashed grey line were studied. If two stars 
were closer than 2" and similar in color and magnitude one of 
the stars was deleted from the sample. In Figure [8] we show 
the resulting contour map. The average background level is 
1.21 • 10~' stars arcmin -2 . We detect overdensities several <x 
above the mean background. These features show sharp edges 
and a flat central plateau. The GC has a pronounced extrati- 
dal halo, but is mostly contained within its Jacobi radius. In 
the northeast we observed distorted contours. Comparable fea- 
tures might be visible on the cluster's southwestern side, but on 
larger scales no features con nected to NGC 2419 are observed. 
ICasetti-Dinescu et al.1 (l2009l) speculated that NGC 2419 might 
be the nucleus of a (former) dwarf galaxy with the Virgo Stellar 
Stream being its tidal tail. 



5.2. NGC 4147 

From the position in the vital diagram (in Figure [1]) we expect to 
observe some hint for the cluster's dissolution due to relaxation. 
NGC 4147 is currently about 18.7 kpc from the Galactic center, 
close to apogalacticon R apo = 25.5 kpc (D99). In Figure [9] we 
plot the contour map of the cluster with a weighted background 
level of n^ k = 3.89 ■ 10 -1 stars arcmin -2 . The large area plot 
reveals a variable background. Overdensities of several cr above 
n^ k are detected. The close-up reveals a complex multiple arm 
morphology. The lcr contour appears to indicate S-shaped tidal 
arms extending over several tidal radii in southern and northern 
direction. We detect a halo of extratidal stars and for the three 
lowest contours we observe multipl e tidal arms. Such a m ulti- 
ple arm morphology is predicted by Montuori et al. (2007J) for 
clusters on elli ptical orbits (e^A j = 0.73) close to apogalacti- 
con. N GC 288 dLeon et al.ll200oi) and Willman I dWillman et al.l 
2006) have a similarly complex morp hology, showing three tidal 
arms. iMartfnez Delgado et al.l (120041) investigated the CMD of 
NGC 4147 and of a comparison field close by. They detected 
some hints for extratidal members stars. In the upper right plot 
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Fig. 5. Number density profiles of NGC 2419, NGC 4147, NGC 5024, NGC 5053, NGC 5272, and NGC 5466. The black open 
squares show the combined profiles. The black filled squares were derived in this study. The light grey solid line is the best fit King 
profile. The dark grey line is the best fit Plummer profile. The horizontal line marks the measured, mean background n^. 



of Figure|5]we plotted the number density profile of NGC 4147. 
A break in the profile's slope is observed. 

iBellazzini et al.l (120031) found stars of the Sagittarius tidal 
arm in the vicinity of NGC 4147. They conclude that this clus- 
ter might be connected to the Sa gittarius dwarf (Sgr dwarf). 
Recently. lLaw & Maiewski (2010) investigated the association 



of MW GCs with the Sgr dwarf based on the authors' newest 
Sgr mode l. The authors conclude, in contrast to earlier stud- 
ies (e.g.. IBellazzini et alJ l200l iForbes & Bridges! 12010). that 
NGC 4147's association with the Sgr dwarf is of relative low 
confidence. It is thus unclear whether the extended features we 
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Fig. 6. Number density profiles of NGC 5904, NGC 6205, NGC 6341, NGC 7006, NGC 7078, and NGC 7089. The black open 
squares show the combined profiles. The black filled squares were derived in this study. The light grey solid line is the best fit King 
profile. The dark grey line is the best fit Plummer profile. The horizontal line marks the measured, mean background ribkg ■ 



observe are a feature of the Sgr dwarf field population. In any 
case, these stars are located at a similar distance as the cluster. 



5.3. NGC 5024 (M 53) 

NGC 5024 is, in projection on the sky, a very close neighbor 
of NGC 5053, although the two clusters are more than 1 kpc 
away from each other in radial direction. NGC 5024 is currently 
located about 18.4 kpc from the Galactic center close to its peri- 
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Fig. 7. Number density profiles of Pal 1, Pal 3, Pal 4, Pal 5, and Pal 14. The black open squares show the combined profiles. The 
black filled squares were derived in this study. The light grey solid line is the best fit King profile. The dark grey line is the best fit 
Plummer profile. The horizontal line marks the measured, mean background n^. 



galacticon R per i = 15.7 kpc (D99). In the middle left plot of 
Figure[5]its number density profile is shown. The King model fit- 
ted to the combined profile reveals extratidal stars, a mild break 
in the slope is visible. In Figure [10] we plot the contour map of 
the cluster. In the left panel, the contour map of the entire ob- 



served area is shown. The cluster is clearly detected as the high- 
est density peak. In the surrounding field mostly features up to 
3cr above the mean of 0.6 stars arcmin~ 2 are observed. Southeast 
of NGC 5024 the close neighbor NGC 5053 is detected, but no 
connection between the two is observed in contrast to the find- 
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Fig. 8. Le/f panel: Contour plot of a large area centered on NGC 2419. Right panel: Blow-up of the contour plot centered on 
NGC 2419. In both panels the isodensity lines are drawn at levels of 1, 2, 3, 5 ,7 ,8 ,10, 20, 40,. . . • cr M? above the mean weighted 
background level. The l<x M( , and 2<Jbkg contours are drawn in light grey, the higher levels are drawn in black. The solid grey arrow 
points towards the Galactic Center. The grey solid line marks the edge of the SDSS survey area. The red solid circle shows the King 
tidal radius. The blue dashed circle shows the Jacobi radius. 
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Fig. 9. Same as Figure[8]but for NGC 4147. The grey dashed arrow indicates the direction of its Galactic motion based on the proper 
motions stated in Table Q] 



ings of IChun et aL ( 201 Oh. The bla ck solid polygon shows the 
area observed by Chun et aT] d2010h . We find a strong (9(f) over- 
density at (198.2, 16.3). In Figure Q~T]this overdensity is not rec- 
ognizable anymore. Therefore, we argue that it is an overdensity 
at the distance of NGC 5024. We found exactly one possible 
BHB/RR Lyrae star centered on the overdensity. It is still de- 
bated in the literature whether NGC 5024 is associated with the 
Sgr ti dal stream (see [Bellazzini et alj|2003bt lLaw & Maiewskil 
2010). Our data do not permit us to conclude whether the over- 



density is related to NGC 5024 or the Sgr stream. We found no 
obvious large scale structure connected with NGC 5024. The 
zoom-in on NGC 5024 shows interesting features. The contours 
are spherical and smooth not only in the cluster center, but also 
in the pronounced extratidal halo. The lcr cont our is strongly 
influenced by the background. IChun et alJ ([2010) also observed 
NGC 5024 and NGC 5053. They detected a tidal-bridge like fea- 
ture around the two GCs and tidal features that may be due to 
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Fig. 10. Same as Figure Ubut for NGC 5024. The large o verdensity in the left panel at (199,17.8) is the GC NGC 5053. The solid 
black line denotes the area observed bv lChun et al. (201(|. 



dynamical interaction between the two clusters. Their findings 
are in contrast to our observations. 



5.4. NGC 5053 



NGC 5053 located about 16.9 kpc from the Galactic center. 
lLauchner etafl (120061) studied the 2d distribution of NGC 5053 
and detected extratidal features east and west of the clus- 
ter. NGC 5053 i s belie ved to be a member of the Sgr dSph 
(Bell azzini et alJ l2003bh . In Figure QT| we show its contour 
map. In the large area contour map we see no pronounced 
large scale structure, but we definitely detect the GC and also 
random background noise up to 3cr above the mean of 4.1 • 
10~ 2 stars arcmin . The large overdensity northeast of the clus- 
ter is NGC 5024. The right panel of Figure QT| is a zoom-in 
on NGC 5053. It reveals distorted contour lines. We observe 
a halo of extratidal stars. Towards the northeast and west the 
contours show a strong asymme try. Similar features were ob- 
served by lLauchner et alJ (120061) . We see a one-armed tail-like 
extension toward the north-west. Such a structure is also seen 
to emerge from NGC 5024 (both in Figures \W\ &GJ}. Our data 
do not permit us t o tell whether the s e curi ous one-armed fea- 
tures are tidal tails. lLaw & Majewskil(l2010h classify NGC 5053 
as moderately likely associated with the Sgr dwarf. Therefore, it 
is possible that we are looking at substructure in the Sgr dwarf 
field population. lForbes~& Bridges (120101) discuss the possibility 
that NGC 5053 or NGC 5024 might be the remnant nucleus of a 
dwarf galaxy. The number density profile is shown in the middle, 
right panel of Figure [5] It is interesting to note that the density 
profile (Figure|5]l do es no t revea l a pronounced extratidal feature. 
Also[Lehmann & S cholzl (Il997l) . o n photographic plat es, did not 
find any indications for tidal tails. IChun et alJ {2010) observed 
NGC 5053 together with NGC 5024, and as stated earlier we do 
not reproduce their findings. 



5.5. NGC 5272 (M 3) 

NGC 5272 is 11.5 kpc away from the Galactic center, located 
close to apogalacticon R apo = 13.7 kpc. Compared to other clus- 
ters in our sample NGC 5272 has a low expected destruction 
rate (G097). The lower left panel of Figure [5] shows its number 
density profile. The combined profile is not well fit with a King 
profile. The profile shows no change in slope as a hint for an 
extratidal halo. In Figure [12] we show the contour map of this 
cluster. Generally, the contours are smooth. On the other hand 
the outer contours up to 3<x show distortions. The extreme dis- 
tortions of the lowest contour are likely to be induced by the 
background. There is no large scale tidal structure observed. The 
background shows random density p eaks of up to 3<x above the 
mean of 4.45 • 10~ 2 stars arcmin -2 . iLeon et all d2000t) studied 
the 2d distribution of NGC 5272 and found extratidal features 
correlated with dust emission. Our contours do not show any 
correlation with the dust. Also Grillm air & Johnso in- 
vestigated the 2d-structure of NGC 5272 with SDSS data and 
did not detect any tidal structure. 

5.6. NGC 5466 

iBelokurov etail d2006l) studied the 2d distribution of NGC 5466 
with S PSS data and detected a 4° -lo ng tid al tail. Earlier stud- 
ies by lOdenkirchen & Grebel d2004l) and iLehmann & Scholzl 
( 1997), on photographic plates, already found indicati ons for the 
tidal tails of this cluster. Grillmair & Johnson (2006) extended 
this tidal tail to 45°. Also. lChun et alJ d2010l) "observed the tidal 
tails of NGC 5466. We also detect the large extratidal feature 
of NGC 5466 (Figure [T3l). The northwest to southeast extension 
of the tidal tail is clearly visible. The contour plot of the larger 
area also reveals the tidal tails close to the cluster. The long ex- 
tensions are hard to distinguish. There are random foreground 
peaks of up to 5<x above the mean of 4.5 ■ 10~ 2 stars arcmin" 2 . 
This is also the case in the plots of lGrillmair & Johnson! J2006). 
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Fig. 11. Same as Figure[8]but for NGC 5053. The large overdensity in the left panel at (198,18) is the GC NGC 5024. 
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Fig. 12. Same as Figure[8]but for NGC 5272. 




Our detections resemble the detections of IChun et al. I d2oToT> . 
NGC 5466 is currently at a distance of 15.8 kpc from the 
Galactic center, close to its perigalacticon R per i = 6.6 kpc (D99). 
Th e inner tidal tails are ali g ned with the clu s ter orb it as predicted 
by iMontuorietaLl (120071) . iFellhauer etaD (120071) modeled the 
destruction of NGC 5466 and were able to reproduce the de- 
tected tidal tails. They further state that this cluster is mainly 
destroyed by tidal stripping at each perigalacticon, and that the 
destruction due to 2-body relaxation plays a minor role. This is 
in agreement with the destruction rates calculated by G097. 



5.7. Palomar 5 

Pal 5 is located 18. 6 kpc from the Galactic center. 
Odenkirc hen et alj d2003l) studied the distribution of cluster stars 
with SDSS data and found tidal ta ils s panning 10° on the 
sky. Alreadv lOdenkirchen etaD (l200ll) an dlRockosi et al.l(l2002h 
found Pal 5's tidal tails with SDSS data. iGrillmair & Dionatosl 
d2006h extended these tails to ~ 20°. lOdenkirchen et alJ d2002h 
measured the line-of-sight velocity dispersion of Pal 5 and con- 
cluded that the dispersion and the surface density profile are 
consistent with a King profile of Wo = 2.9 and r, = 16.1', 
similar to t he values found here . The disruption of Pal 5 was 
modele d in iDe hnen et al. (2004) base d on the initial observa- 
tions of lOdenkirchen et all 02003). Viv as & Zinnl d2006h studied 
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RR Lyrae stars in the QUEST survey. They observed an overden- 
sity of RR Lyrae stars in the vicinity of Pal 5. Some of these stars 
trace the cluster's tidal tails. The middle right panel of Figure [7] 
shows the number density profile of Pal 5. The overabundance 
of member stars due to the tidal tails is clearly visible as a pro- 
nounced break in the profile. In Figure [14] we show the derived 
contour maps. In the large area view the tidal tails are clearly 
visible emanating from the cluster. In the zoom-in the «S-shape 
is observed as predicted by theory. The background harbors fluc- 
tuations several cr above the mean of 1.1 • 10 _1 stars arcmkr 2 . 
iKoch et ail (120041) show that the mass segregation in Pal 5 ex- 
tends to its tidal tails, which contain predominantly low mass 
stars. Based on a kinema tic study of the stars in the tails of Pal 5, 
lOdenkirchen et al.1 (2009) suggest that the cluster's orbit is not 
exactly aligned with its tails. 

5.8. NGC 5904 (M 5) 

NGC 5904 is currently 6.1 kpc away from the Galactic center, 
close to its perigalacticon R pe ri = 2.5 kpc (D99). I.e., poten- 
tial tidal tails would be expected to be aligned with the clus- 
ter's orbit. In the number density profile, shown in Figure [6j 
no change in slope is observed. In Figure [TBI we show the con- 
tour map. The field background shows contours at levels of 2<x 
above the mean of 2.1 • 10 _1 stars arcmkr 2 . The background 
density increases towards the Galactic center, but a connection 
to NGC 5904 does not seem to exist. The tidal tails of Pal 5 
are not directly visible, only the main body of this neighboring 
cluster is detected as an overdensity at (229, 0). The three lowest 
contours of NGC 5904 are distorted, the lcr-contouris compara- 
ble to the background, while the 2,3,4<x-contours indicate an ex- 
tratidal halo. The contours are elongated in the directi on of mo- 
tion p ossibly indicating weak evidence of tidal tails. iLeon et ail 
(2000) also studied the 2d distribution of this cluster. They also 
observed round contours, but no extratidal halo. NGC 5904 is 
appr oaching perigalacticon, s tarting its crossing of the Galactic 
disk (Odenkir chen et al. Ill997t) . 



5.9. NGC 6205 (M 13) 

NGC 6205 is currently at a distance of 8.2 kpc from the Galactic 
center, close to its perigalacticon R per i - 5.0 kpc (D99). Thus, 
we expect potential tidal tails to be aligned with the orbital path. 
In Figure[16]we show the resulting contour plot. The background 
shows residuals up to 2<x above the mean of 1.2 ■ 10~' stars 
arcmkr 2 . No large scale features are observed. The cluster's 
central contours are smooth. The contours in the outer regions of 
the cluster show distortions. We detect a halo of extratidal stars 
with a slight elongation in the direction of motion. The num- 
ber density profile in Figure [6]revea ls an o verabundance of stars 
around the tidal radius. ILeon et I] §000) studied the 2d distri- 
bution of NGC 6205 as well. They also detected disturbed con- 
tours in the southeast, corresponding to the feature visible in our 
data. 



5.10. NGC 6341 (M 92) 

NGC 6341 is located 8.2 kpc away from the Sun, close to its 
apogalacticon, R apo - 9.9 kpc (D99). In Figure[T7]we show our 
derived contours of NGC 6341. In the left panel the effect of 
missing data due to the boundary of the SDSS scan regions in 
the southeast corner is visible. Northeast of the cluster we de- 
tected an artifically doubled background density. We corrected 
for this. All stars in the area denoted by the dashed grey line were 
studied. If two stars were closer than 2" and similar in color and 
magnitude one of the stars was deleted from the sample. As seen 
in the contour map in this way we were able to reach almost 
the same background level as in any other area of our contour 
map. The mean background level was determined in the area 
northeast and southwest of the cluster ignoring the spurious ar- 
eas. The mean field background is 9.6 • 10~ 2 stars arcmkr 2 . The 
background shows overdensities with contours of 3<x and more 
mainly southwest of the GC. 

The zoom-in in the lower panel of the same figure reveals 
mostly smooth contours in the cluster center. The main distortion 



14 



Katrin Jordi and Eva K. Grebel: Search for Extratidal Features Around 17 Globular Clusters in the Sloan Digital Sky Survey 




is in the direction of the corrected background and due to the 
missing data on the opposite side of the cluster it is not possible 
to say whether this feature indicates a tidal tail. But also along 
a southwest to northeast axis we detect extratidal m aterial. The 
2d stru cture of NGC 6341 has been studied before by Test a et ail 
( 2000). They observed somewhat elongated contours in the same 
direction as we do. The number density profile in Figure[6]shows 
an overabundance of stars around the tidal radius. 

5.11. NGC 7006 

NGC 7006 is one of the most remote clusters in our sample with 
a distance to the Galactic center of Rmw = 38.8 kpc. The cluster 
is currently moving toward the Galactic plane and the Galactic 



center dDinescu et al. I l200lh . It has not been included in any 
study on tidal tails searches so far. In Figure[18]we show its con- 
tour map. We clearly detect the GC as the highest overdensity. 
The lcr-contour is very distorted and comparable to the back- 
ground noise. The mean field background level is 1.2 ■ 10~' stars 
arcmin~ 2 . The 2<x and 3<x contours show some symmetric dis- 
tortion. All higher contours are smooth, although the halo of ex- 
tratidal stars for NGC 7006 is huge. For no other cluster in our 
sample we find such a large extratidal halo. 

5.12. NGC 7078 (M15) 

NGC 7078 is currently located 9.8 kpc from the Galactic cen- 
ter, close to its apogalacticon R apo = 10.4 kpc (D99), i.e., we do 
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not expect any correlation between the orbit and the (possible) 
tail's orientation. In Figure [19] we show the contour map. The 
field around the cluster is very smooth. It shows two peaks of 3<x 
above the mean of 3.0 • 10 _1 stars arcmkr 2 , but not in any sym- 
metric configuration around the cluster. No large scale structures 
are detected. The contours at the cluster center are undisturbed, 
but the contours at the tidal radius show some distortions. The 
proximity of NGC 7078 to the edge of the SDSS survey intro- 
duces some asymmetry. Neverthele ss, we observe an extr atidal 
extension in southwestern direction. Grillm air et alj (fl995) stud- 
ied the 2d distribution of NGC 7078. They found an excess of 
cluster stars extending toward the southeast and a further density 



peak just north of the cluster. These features do not show up in 
our data. Instead, we confirm the features seen in the immediate 
surro undings of NGC 7078 that were also found bv lChun et al.l 
(2010). Owing to the borders of the SDSS survey, their north- 
western feature is not co vered by our dat a. On a larger scale, out- 
side the field of view o flChun et alJ d2010h . we see extended low- 
density (1 - 2<t) structures in roughly east-west direction. These 
are not aligned with the direction toward the Galactic center, but 
appear at an angle of approximately 20 - 30 deg (Figure [T9l , 
suspiciously aligned with the direction of the SDSS scan. While 
suggestive, deeper data are needed to confirm or rule out whether 
these are extended tidal tails. We observe a similar extratidal 
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Fig. 17. Same as Figure[8]but for NGC 6341 . The solid grey lines denote the edge of the SDSS footprint. The dashed grey line marks 
the area where the stellar density was doubled (see text for details). 



structure of the contourlines. Figure|6]shows the number density 
profile of NGC 7078. We detect a change in the profile's slope. 

5.13. NGC 7089 (M2) 

NGC 7089 is located 10.1 kpc from the Galactic center, close to 
perigalacticon R per j = 6.4 kpc (D99). I.e. a possible inner tidal 
tail is expected to be aligned with the orbital path of the cluster. 
In Figure [20] we show the contour plot of NGC 7089. The inner 
contours are smooth. The extended lcr-contouris comparable to 
the background contours. The overall background shows at most 
contours of 2<x. The mean background level is at 1.2 • 10 _1 stars 
arcmin~ 2 . No large scale features are detected. NGC 7089 shows 
an extratidal halo, which is a sp herical structure. No tS-shape 
or elongation is observed. Also iGrillmair et al.l (1 19951) studied 
the 2d structure of NGC 7089. They found an excess of cluster 
stars extending along an E-W axis. In our contours only the 2cr- 
contours might show some elongation, but overall the halo is 
quite circular. 

iDalessandro et al.l ([2009) derived a number density profile 
and also detected a slight excess of stars beyond the tidal radius 
compared to the best fit King model. Our number density profile 
is shown in the Figure|6] We detect a change in slope around the 
tidal radius as expected for a halo of extratidal stars. 

5.14. Palomar 1 

Pal 1 is not a typical GC. Its CMD is quite unusual. It 
has a populated main sequence, but hardly any stars on the 
re d giant and the horizon t al bra n ch. This feature was no ticed 
bv iBorissova & Spassoval d 19951) . iRosenberg et alJ dl998h . and 
ISaraiedini et al.l d2007l) . Further, this clus ter might be up to 8 Gy r 
younger(!) than a typical GC like 47 Tuc (ISaraiedini et al.l2007b . 
It is also discussed in the literature, whether Pal 1 is a member 
of the Monoceros stream /Canis Major dwarf (ICrane et al . 2003; 
Forbes & Bridges 20Tob . It cannot be ruled out that Pal 1 is mis- 
classified as a GC since objects of this young an age are usually 



considered open clusters. The low background density contour 
map shown in Figure [21] is due to the cluster's high declina- 
tion plotted in Galactic longitude and latitude, (I, b). Although 
Pal 1 has a high destruction probability (G097), mainly driven 
by disk and bulge shocks, it does not show any distortions of its 
contours, except for a large circular halo of extratidal stars (com- 
parable to NGC 7006). In the field we see no density peaks. The 
number density profile is shown in Figure [7] The profile drops 
outside of ~ 30 arcmin because the edge of the survey is reached. 
In summary we can say that we detect a symmetric overdensity 
at the position of Pal 1 . 

5.15. Palomar 3 

Pal 3 is currently at a distance of Rmw = 84.9 kpc, very close 
to perigalacticon R per i - 82.5 kpc (D99). As this is the clus- 
ter's closest approach to the Galaxy we do not expect to see any 
large signs of dissolution due to tidal shocks. Also the destruc- 
tion rates derived by G097 are small. In Figure[22]we show the 
contour map of a large area centered on Pal 3. The cluster can 
clearly be detected as the strongest overdensity. Nevertheless, 
the background shows several pronounced overdensities. The 
largest overdensity at (153.5,-1.5) is the signal of the dwarf 
spheroidal galaxy Sextans. The cluster's contours look undis- 
turbed. We observe a distortion of the contours towards the 
Galactic center, but no corresponding feature on the opposite 
side. Generally, Pal 3 has a large extratidal halo which shows 
a sharp edge. The number density profile is shown in Figure [7] 
We note that owing to the sparseness of Pal 3's CMD and its 
large distance deeper data are needed to investigate options such 
as a connection with Sextans. 

5.16. Palomar 4 

A density plot of Pal 4 and surroundings is shown in Figure |231 
We discover a large halo of extratidal stars. The lcr-contour is 
stretching along towards the Galactic center, but only in one di- 
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Fig. 18. Same as Figure[8]but for NGC 7006. The solid grey lines denote the edge of the SDSS footprint. 



rection. The field around Pal 4 is showing strong overdensitie of 
up to lOcr. These density peaks show steep edges and in the cen- 
ter broad plateaus. It is unclear where these features are coming 
from. Comparing these field peaks with the cluster's tidal tail, 
hardly any difference can be seen. An interesting structure of the 
contour lines is visible at the cluster's western edge. 

ISohn et al.l d2003l) studied the 2d-morphology of Pal 3 and 
Pal 4 using CFH12K wide-field photometry. For Pal 3 they de- 
tect some extensions in the direction towards the Galactic center 
and anti-center which are not more prominent than lcr over the 
background. For Pal 4 they see a tail extending to the NW. In our 
data we do not observe such a feature. Pal 4 is the most distant 
cluster in our sample. Again, deeper data are needed to reveal 
the possible existence of tidal tails. 



5.17. Palomar 14 



In Figure l24l we show t he contour map of Pal 14. 
iMartmez Deleado et al.l (l2004h studied the CMD of Pal 14 and 
compared it to the CMD of a region outside the cluster. They de- 
tected some hints for extratidal stars in a region of the CMD. 
We see smooth central contours and an extratidal halo. The 
field background shows some overdensities of up to 5cr above 
the mean which is 1.47 • 10 _1 stars arcmin~ 2 . The lcr and 2<x 
contours show an elongated distortion, comparable to the field 
background contours. The extratidal halo is spherical with some 
asymmetry towards the southeast. Whether this indicates the be- 
ginning of sparsely populated tidal tails remains unclear. Also 
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Fig. 19. Same as Figure[8]but for NGC 7078. The solid grey lines denote the edge of the SDSS footprint. 



for this GC, deeper data would be desirable in order to confirm 
this or to rule it out. 



6. Discussion 



6. 1 . Horizontal branch stars 
From theoretical investigations 



g al 

iBaumgardt&Makinol H003) we 
stars of lower mass are lost. Measurements of the mass functio n 
of Pal 5 and its tidal tails confirm this picture dKoch et al.l l2004). 
Also the number density of main-sequence stars is larger than 
for RGB stars. Therefore, tracing tidal structures with only 
upper RGB stars may be difficult, as is the case for Pal 3, Pal 4, 
Pal 14 and NGC 2419 and this es exacerbated by the sparsely 



and 
know 



simulations (e.g., 
that preferentially 



populated RGBs of the Palomar GCs. Nonetheless, we do find 
halos of extratidal stars as possible signs for some dissolution 
for all four clusters. But the field background is noisy and shows 
overdensities of more than 3cr. 

The blue horizontal branch (BHB) and RR Lyrae stars are 
located at a position in the CMD where hardly any field stars are 
located. Therefore any BHB or RR Lyrae star at the cluster dis- 
tance found outside the cluster's tidal radius is wi th a high proba- 
bility a former cluster member. For NGC 5466, Belokurov et al.l 
(2006) found that the BHB/RR Lyrae stars also trace the tidal 
tails. For those of our GCs that show no large tidal structures, but 
only distorted tidal contours the majority of the BHB/RR Lyrae 
stars are inside the tidal radius. As an example: In NGC 6341, 
we found 194 BHB/RR Lyrae stars. Only 20 of those are found 
outside the tidal radius. The remaining 174 BHB/RR Lyrae stars 
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Fig. 20. Same as Figure[8]but for NGC 7089. The solid grey lines denote the edge of the SDSS footprint. 



are all found within the tidal radius (red ellipse in Figure ITTb . 
For NGC 4147, for which we have found a complex multiple 
arm morphology and extended lcr structures that look like weak 
tidal tails, most BHB stars are within the innermost drawn con- 
tour, and 30% are found outside. The BHB stars do not trace 
the multiple arms. But the BHB/RR Lyrae stars in the field are 
mostly found in overdensities. E.g., at (181.2, 17) two overden- 
sities of 9<x and 10cr, respectively are found. In the center of 
both a single BHB/RR Lyrae star is detected. At the same time, 
there are BHB stars that are not connected to any overdensity. 
The algorithm used to derive the density of cluster member stars 
gives high weights to BHB and RR Lyrae stars. Without spec- 
troscopic information it is not possible to draw any reliable con- 
clusion whether BHB/RR Lyrae stars connected to a prominent 



overdensity in the background are 'lost' cluster members. But 
the existence of BHB/RR Lyrae stars without any overdensity is 
supporting the hypothesis that the stars in those overdensities are 
more likely former member stars. 

A similar connection between some of the cluster's 
BHB/RR Lyrae stars and overdensities in the field is observed 
for NGC 5024, NGC 5053, NGC 5466, NGC5904, NGC 6205, 
and NGC 7006. Our current data do not permit us to conclude 
whether these overdensities are related to the GCs. 

For some GCs in our sample a connection to the Sagittarius 
stream (e.g., lBellazzini^t~aTI l2003bt lLaw & Maiewskil 1201 



Forbes & Bridges 2010) or the Monoceros stream (Crane et al 
2003: [Forbes & Bridgesll2010l) is discussed in the literature. It is 
possible that the detected extratidal structures and overdensities 
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in the field of NGC 4147, NGC 5024, NGC 5053 and Pal 1 are 
features induced by the field populations of the disrupted dwarf 
galaxies. The data used here do not allow us to make a definite 
conclusion. 



6.2. Comparison to theory 

The vital diagram of G097 in Figure Q] predicts a dissolution of 
NGC 4147 due to 2-body relaxation. The multiple arm feature 
observed in Figure [9] for NGC 4147 sup ports the cluste r 's dis- 
solution. Furthermore, the simulations bv lMontuori et al.l (120071) 
found such a multiple arm morphology for GCs on eccentric or- 
bits approaching apocenter. NGC 4147's current position close 
to apogalacticon matches the simulations well. Of all stars in the 



distorted lcr-contour around NGC 4147 about 13% are found 
outside the 4cr-contour. The 4cr-c ontour corresponds roug hly to 
the Jacobi radius of NGC 4147 dBaumgardt et alj 12009). i.e., 
stars outside the 4<x-contour left the potential field of the cluster. 
The tidal features of NGC 4147 are small compared to the ones 
for Pal 5. Pal 5 is close to destruction due to tidal shocks, while 
NGC 4147's dissolution is mainly driven by relaxation. 

The Palomar clusters in our sample also lie outside the tri- 
angles in Figure Q] As the outermost triangle is only appropriate 
for clusters with a current distance to the Galactic center smaller 
than 12 kpc, and Pal 3, Pal 4, and Pal 14 are currently further 
than 69 kpc away from the Milky Way no conclusions can be 
drawn from the diagram for these three clusters. Pal 5 is known 
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Fig. 23. Same as Figure[8]but for Pal 4. 



f5 



I i 



Pi ° 



Vij 




1 t 
. . 



do) '-c 




Fig. 24. Same as Figure[8]but for Pal 14. 



to have tidal t ails. Its destruction rate was studied in detail by to detect tidal tails if present. The question of whether Pal 1 is 
iDehnen et alj (120041) . really a GC or rather an open cluster remains unsolved. 



Pal 1 is also located outside the triangles in the vital dia- 
gram. Hence a dissolution due to 2-body relaxation is likely, but 
the contour map in Figure I2TI shows no typical .S-shape. It re- 
veals a large, spherical extratidal halo though. The halo extends 
much further out than the cluster's Jacobi radius. The Jacobi ra- 
dius is slightly smaller than the 20<x-contour (the second most 
inner contour). 71% of all stars located within the lcr-contour 
are found outside the Jacobi radius, i.e., should not be bound 
to the GC anymore. Pal 1 in general is an interesting object. Its 
CMD shows a clear main-sequence but a very sparsely popu- 
lated red giant branch. The low stellar density makes it difficult 



The remaining GCs in our sample from the NGC catalog 
are all located well within the triangles of G097's vital di- 
agram. Hence, no (pronounced) signs for dissolution are ex- 
pected. NGC 5466 is an exception and shows that the theoretical 
calculations can be improved, as tidal tails have been detected. 
For NGC 5053 some contours are located outside the Jacobi ra- 
dius, a possible sign for unbound cluster stars in the northeast, 
west, and southwest. Generally the clusters show an extratidal 
halo, which is to be expected from dynamical studies. These 
GCs are not in danger to be destroyed within the next Hubble 
time, according to G097 and D99. This is supported by the lack 
of large scale tidal features in our study. 
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Table 4. Power law slopes Q b a of the radial density profiles be- 
tween the two radii a ■ r, and b ■ r t . The stated errors are only the 
errors from the fit. 



NGC Name 



-G? 



-Gf 



2419 




3.44 ± 0.75 


3.32 


±0.47 


1.13d 


:2.4 


4147 




1.48 ±0.24 


1.70 


±0.24 


4.19 ± 


0.32 


5024 


M53 


3.73 ±0.35 


3.69 


±0.24 


1.68 ± 


1.28 


5053 




0.33 ± 1.05 


2.35 


± 1.76 






5272 


M 3 


0.94 ± 0.38 


0.91 


± 0.16 






5466 




3.89 ± 0.46 


1.52 


±4.01 


3.39 ± 


0.75 


5904 


M5 


3.33 ±3.72 


3.33 


±0.24 


-0.70 i 


= 1.83 


6205 


M 13 


3.28 ± 0.33 


3.28 


±0.20 


1.65 ± 


6.80 


6341 


M92 


3.65 ±0.21 


3.76 


±0.21 


3.76 ± 


11.0 


7006 




3.38 ±0.73 


3.36 


±0.51 


1.23 ± 


6.38 


7078 


M 15 


2.88 ± 0.26 


3.00 


±0.21 






7089 


M2 


2.58 ±0.21 


2.65 


±0.15 


5.09 ± 


2.58 




Pal 3 


7.14 ± 10.58 


7.16 


±8.27 








Pal 4 


1.47 ± 1.15 


1.88 


±0.77 








Pal 5 


0.26 ± 0.42 


1.33 


±0.29 


3.48 ± 


2.95 




Pal 14 


3.13 ±2.96 


3.46 


±2.58 







6.3. Number density profiles and halos of extratidal stars 

Most of the clusters in our sample show a (pronounced ) halo 
of stars outside the tidal boundary. ICombes et al.l d 19991) sim- 
ulated GCs in a tidal field and observed such a halo as well. 
The stars outside the tidal radi us spread out in density like a 
power law. ICombes et ail (1 19991) transformed the 3 -dimensional 
simulations into a (observable) projection on the sky and de- 
rived the power law to be "r~ 3 or stee per". The extratidal ha- 
los observed by iGrillmair et all d 19951) a nd iLeon et all (2000) 
show only slopes of -1. Combes et alj d 1999b argue that the 
discrepancy is a consequence of noisy background-foreground 
subtraction. Our stellar sample contains very little contamina- 
tion. Therefor e, an investiga t ion o f these halos is of interest. 
Furthermore, Grill mair et all (1 19951) observed for clusters with 
tidal extensions a break in their number density profiles, becom- 
ing pure power laws a t larger radii. 

ILeon et ail (2000) introduced Q h a as the slope of the number 
density profile between a ■ r, and b ■ r,. They chose to compute 
the three slopes: Q^, Q h v and Q 6 y . We applied the same method 
to the clusters in our sample and derived the three slopes. In the 
profile shown in Figure [5]|7] slope changes are visible in many 
clusters. The slope was only computable for clusters, where 



the background was not too noisy. In order to fit the power law 
slopes we performed a least-squares fit between the two radii. In 
Table|4]we list the resulting slopes for all clusters. 

The theoretical prediction of slopes of "r~ 3 or steeper" is ob- 
served for many clusters in our sample for . Only NGC 4147, 
NGC 5053, Pal 4, and Pal 5 have a significantly flatter slope than 
predicted. The measurement for Pal 3 is affected by a large error. 
The fits for show large errors or were in some cases not pos- 
sible at all, due to the noisy background. For Pal 1, we did not 
perform any fit due to the strange shape of the number density 
profile. In general, the simulated halos of ICombes et all d 19991) 
reproduce well the observed halos. 

For most GCs an extratidal halo exists, while the Jacobi ra- 
dius is larger than this halo. This is the case for NGC 2419, 
NGC 5024, NGC 5272, NGC 5904, NGC 6205, NGC 6341, 
NGC 7006, NGC 7078, and NGC 7089. From the definition of 
the Jacobi radius this suggest that the stars in the extratidal halo 
have not left the cluster potential. The derivation of the Jacobi ra- 
dius is based on measurements and assumptions, e.g., the mass- 



to-light ratio is assumed to be 2, which may not be generally 
true, see M05 for mass-to-light ratios of the Galactic GCs. One 
possibility is that the contour maps do not reveal halos of extrati- 
dal stars_bir^jti^theJialo of the GCs. Although, the simulations 
bv lCombesetal] (fl999) reveal extratidal halos with the identical 
power-law slope as our observed halos. 

The Palom ar clusters have a h alo of stars outside the 
Jacobi radius. iBaumgardt et al.l d2009l) classifies these clusters as 
"tidally filled". We actually observe that these clusters are even 
more extended. Pal 14 is a very special case: the Jacobi radius 
is smaller than the tidal radius from the King profile. The low 
measured v elocity dispersio n of Pal 14 does not suggest a dis- 
solved GC dJordi et alJ2009l) . neither do the mostly symmetrical 
contour lines. For Pal 5 the same observation is made, the Jacobi 
radius is smaller than the tidal radius. Pal 5 has a lower than ex- 
pected velocity dis persion, although it will b e destroyed within 
the next 110 Myr dOdenkirchen et al.|[2002h . It is possible that 
the large distance to Pal 14 prevents us from detecting its tidal 
tails although we do see possible extensions at the 1 - 2<x level. 
Or their orientation makes them difficult to detect. Simulations 
of dwarf spheroidals have shown that the detection of the char- 
acteristic 5-shape of tidal tails depends on th e position of the 
orbit relative to the line of sight (Munoz et al. 20081). 

7. Summary & Conclusions 

We studied the 2d-distribution of (potential) member stars on 
the sky of the 17 GCs in the SDSS. The stars were chosen from 
the SDSS DR 7 and from An08. All stars in a field around a 
GC were selected in color and magnitude to match the observed 
colors and magnitudes of cluster stars. 

We derived the number density profiles for the 17 clusters in 
our sample. Thanks to the large area coverage of the SDSS, we 
were able to trace the profiles as far out as to the flat field star 
background. In some cases a break in the profile was detected 
and interpreted as a sign for an extratidal halo. Comparing the 
profiles with the 2d-distribution of cluster stars on the sky, we 
detected a halo of apparent extratidal stars for more clusters than 
a clear break in the profile. 

A color-magn i tude weighted counting algorithm from 
lOdenkirchen et all d2003l) was used to derive density maps of 
cluster stars on the sky. We recovered the known extended tidal 
tails of Pal 5 and NGC 5466. For NGC 4147 we detected a com- 
plex, two arm feature of extratidal stars, and for NGC 5904 we 
find a roughly symmetric extension that might indicate short 
and stubby tails. For NGC 7006 and Pal 1 the extratidal halo 
is remarkably large. For Pal 1, 71% of the stars of the GCs are 
found outside the Jacobi radius, i.e., are not bound to the cluster 
anymore. Nevertheless, they are arranged spherically around the 
cluster. No sign of tidal effects is observed. We also detected 
ex tratidal features for N GC 5053, comparable to the findings 
of lLauchner et al.l (|2006). We confirm e xtratidal features i n the 
immediate surroundings of NGC 7078 dChun et alj|2010t) . For 
Pal 14 deeper data are needed in order to investigate whether the 
weakly visible extensions constitute tidal tails. 

The comparison of our observations to the theoretical study 
of the GCs' destruction rates by G097 shows generally good 
agreement. Most NGC clusters in our sample were not expected 
to show signs for dissolution. They mostly only revealed a halo 
of extratidal stars well within their Jacobi radius. These extrati- 
dal halos show for the majority of our sample the theoretically 
expected "r~ 3 " density slope. 
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